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In this talk*, we describe recent experimental progress in detecting the chiral anomaly in the Dirac 
semimetal NasBi in the presence of a magnetic field. The chiral anomaly, which plays a fundamental 
role in chiral gauge theories, was predicted to be observable in crystals by Nielsen and Ninomiya 
in 1983 1 . Theoretical progress in identifying and investigating Dirac and Weyl semimetals has 
revived strong interest in this issue HE]. In the Dirac semimetal, the breaking of time-reversal 
symmetry by a magnetic field B splits each Dirac node into two chiral Weyl nodes. If an electric 
field E is applied || B, charge is predicted to flow between the Weyl nodes. We report the observation 
in the Dirac semimetal NasBi of a novel, negative and highly anisotropic magnetoresistance (MR). 

We show that the enhanced conductivity has the form of a narrowly defined plume that can be 
steered by the applied B. The novel MR is acutely sensitive to deviations of B from E, a feature 
incompatible with conventional transport. The locking of the current plume to B appears to be a 
defining signature of the chiral anomaly. 


I. PREAMBLE 

The intriguing possibility of observing the chiral (or 
axial) anomaly as a novel charge current in a crystal has 
been discussed since 1983 [L. In the field of topolog¬ 
ical phases of matter, this question has lately received 
intense interest in the context of Weyl semimetals [2Hz]. 
To realize a Weyl semimetal, a promising path appears 
to be to start with a 3D Dirac semimetal puns], and 
then split each Dirac node into a pair of Weyl nodes 
using a magnetic field B. Rapid progress in identify¬ 
ing candidates for Dirac semimetals has been made m 
Q2]. The two semimetals identified mm, NasBi and 
CdaAs2, were recently confirmed to be Dirac semimet¬ 
als by angle-resolved photoemission (ARPES) [191123] . 
scanning tunneling microscopy [24] and transport exper¬ 
iments [25] [26] • We focus here on the results from a 
recent experiment in which the axial current in NasBi 
is detected as a collimated current plume locked to the 
direction of b EZ1. Figure [l] summarizes pictorially our 
main finding. 


II. NODE PROTECTION BY SYMMETRY AND 
WEYL NODES 

The problem of how the energy gap in a narrow-gap 
semiconductor closes and re-opens as a function of a 
tuning parameter V (e.g. pressure, doping or temper¬ 
ature) has received considerable attention in the con¬ 
text of topological phases and 3D Dirac states in the 
bulk [T2I1T8] . The situation in which the gap closes at 
a single, critical value V c used to be regarded as the 
norm (Fig. [2^i, b). This situation - dubbed acciden- 


* Article based on a talk by NPO at the APS March Meeting, San 
Antonio, 2015. 



FIG. 1: A narrow current plume steered by magnetic field 

B. 


tal band crossing - is inherently unstable and not very 
interesting experimentally. Recent progress has shown 
that Dirac nodes can be protected by symmetry over 
an extended range of V (and hence stable) (Fig. HD- 
If time-reversal symmetry (TRS) and inversion symme¬ 
try (IS) are both present, the protection extends only 
to Dirac nodes pinned to high-symmetry points (either 
T or the corners of the Brillouin Zone) [14]. Unfortu¬ 
nately, this severely limits the number of realistic can¬ 
didate materials. Later, it was realized that inclusion 
of point-group symmetry (PGS) extends the protection 
to nodes away from the high-symmetry points provided 
they lie on a symmetry axis pm]. This crucial insight, 
which greatly expands the list of materials, led to the pre¬ 
diction by the IOP group [Tbl 17] that NasBi (with the 
PGS C3) and CdsAs2 (C4) are Dirac semimetals with 
protected nodes (Fig. [2J3) . In the past year, ARPES \19f- 
[23;, STM [24] and transport experiments [2_5] [26] have 
confirmed these predictions. Recently, surface FS arcs 
have been detected by ARPES [29] . 

When a Dirac semimetal is immersed in a strong mag¬ 
netic field B, one of the 3 protecting symmetries, TRS, 
is removed. Instead of reopening the gap, each of the 
Dirac nodes is predicted to split into two Weyl nodes H- 
[6] (Fig. [3jA) . The Weyl nodes, which may be regarded 
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FIG. 2: Comparison of band closing scenarios as a function 
of the tuning parameter V (here written as m). The upper 
panel (a) shows the accidental band-closing case in which TRS 
and IS stabilize the Dirac node at a high-symmetry point 
at one critical value m c . The lower panel shows the case 
when TRS, IS and C n protect the Dirac node on a symmetry 
axis over the finite interval m c 1 —)> m c 2 (reproduced from 
Yang and Nagaosa [18]). The two predicted nsnui Dirac 
semimetals NasBi and CdsAs 2 correspond to C 3 and C 4 point 
group symmetries. 


FIG. 3: Panel A: Breaking of time reversal symmetry (TRS) 
in an applied B splits the Dirac node into a pair of Weyl nodes 
with opposite x- The Weyl node with x = 1 (-1) is a monopole 
source (sink) of Berry curvature J-. The Chern flux captured 
by the Gauss surface S (pink sphere) normalized by 2i r gives 
the value of x- Panel B: The Adler-Bell-Jackiw anomaly de¬ 
scribing the very rapid decay of 7r° into two photons 7 m [ 9 ]. 
Panel C: The proposal of Nielsen and Ninomiya pQ for observ¬ 
ing the chiral anomaly in a ID crystal. The slanted lines rep¬ 
resent the dispersion of right- and left-moving chiral fermions 
described by the massless Dirac equation 27 ^ <9^ T = 0. Chi¬ 
ral charge conservation is broken when an E-field is applied. 
Charge is pumped between the two branches. The dashed 
curve indicates the connection between the chiral branches 
deep in the Fermi sea in the solid-state context. 


as monopole sources and sinks of Berry curvature T in 
k space, come in pairs with opposite chiralities x = Ti¬ 
lt is helpful to regard E(k) as an effective magnetic field 
that lives in fc-space. If we surround a Weyl node with a 
Gauss surface <S, the Chern flux captured defines the chi¬ 
rality x = (1/ 27r) ^E(k) • dS( k). Thus, breaking TRS 
produces an intense Berry curvature T in the vicinity 
of each Weyl node. Given that T imparts a transverse 
“anomalous” velocity = J x eE to a wave packet, we 
should expect a host of novel transport features (E is an 
applied E-field). The most interesting involves the chiral 
anomaly. 


III. CHIRAL ANOMALY 

A bit of topological physics fell into quantum field 
theory (QFT) in the late 1960s [SlfTT]. The charged 
pions 7r^~ are remarkably long-lived mesons (lifetimes 
r ^ 2.3 x 10 -8 s) because, being the lightest hadron, 
they can only decay by the weak interaction into muons 
and neutrinos via the processes 7r - fi~ + and 

7r + /i + + Mysteriously, the neutral pion 7r° de¬ 

cays much more quickly (by a factor of 300 million) even 
though it is a member of the same isospin triplet. Instead 
of the slow leptonic channels, 7 r° decays by coupling to 
the electromagnetic field F^ v in the process 7r° 27 . 

The relevant diagram (Fig. [3(3), called the Adler-Bell- 
Jackiw anomaly 0E], is a triangular fermion loop that 
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links the i r° (the axial current) to the 2 photons (vec¬ 
tor currents) [lOj [lT] . A hint of the topological nature 
of the anomaly is that the one-loop diagram receives no 
further corrections to all orders of perturbation theory. 
Subsequent research revealed that the anomaly expresses 
the breaking of a classical symmetry by quantum fluctu¬ 
ations. In modern QFT, the anomaly plays the funda¬ 
mental role of killing unviable gauge theories mm- a 
proposed chiral gauge theory must be anomaly free. Oth¬ 
erwise it is not renormalizable. Arguably the most im¬ 
portant example of the anomaly-free rule is the Glashow- 
Salam-Weinberg electroweak theory, in which the 4 trian¬ 
gle anomalies linking the lepton and quark doublets with 
gauge bosons sum exactly to zero within each genera¬ 
tion. This fortuitous cancellation has been called “mag¬ 
ical” [TO] . 

In 1983, Nielsen and Ninomiya (NN) [1] proposed 
that the chiral anomaly may be observable in a crystal 
with massless Dirac states in even spacetime dimension 
(1+1 or 3+1) (Fig. IT)- When we set the mass m— 
0 in the Dirac equation, the Dirac 4-spinor T describes 
two independent populations of chiral particles with spin 
locked either parallel to their momentum (and described 
by the Weyl 2-spinor i/jr) or antiparallel (V’l)- The chi¬ 
ral charges are separately conserved. This chiral symme¬ 
try is ruined as soon as we couple the particles to 
(by applying an E field). Chiral charge is pumped from 
one branch to the other at a rate given by the triangle 
anomaly, viz. 


W = 


X 47r 2 h 2 


E B. 


(1) 


At first blush, this would appear to be a restatement 
of charge transport in an ordinary metal. If we regard 
the left and right-moving states as opposite limbs of a 
large Fermi surface (FS), the E field simply biases the 
occupation of the two limbs leading to the familiar FS 
shift in k- space. However, as discussed below, the implied 
“locking” of E to B is unusual and new. It actually 
leads to features that fall beyond the purview of standard 
transport in metals. 

Unsurprisingly, the discussion in 2011 of Weyl nodes in 
the context of 3D topological phases of matter EH3 stim¬ 
ulated intense renewed interest in NN’s proposal, initially 
among theorists but now extending to experimentalists 
as well. Theoretically, the states in the two Weyl nodes 
are quantized into Landau levels in strong B. A key fea¬ 
ture of the Weyl spectrum is that the n = 0 Landau level 
(LL) is chiral, with velocity v parallel or antiparallel to B 
(Fig. §\). We have precisely the situation described by 
NN, except that the ID dispersion direction is controlled 
by B 

As mentioned, the breaking of time-reversal symme¬ 
try (TRS) by B causes each Dirac node to split into two 
Weyl nodes [HHH]- In the absence of an electric field 
E, the two chiral populations are separately conserved 
(Fig. §\). However, application of E ruins the conser¬ 
vation and leads to a charge-pumping process between 




20 30 40 50 60 70 80 90 


2 theta (degree) 


FIG. 4: Crystal growth and structure analysis, (a) Photo¬ 
graph of stoichiometric cooled NasBi (1) crystals with nat¬ 
urally cleaved pristine faces, (b) Photograph of typical 90% 
Na-flux grown NasBi (2) crystals, (c) The honeycomb crystal 
structure of NasBi. The Na(2) atoms are cross linked to Bi 
of -Na(l)-Bi-Na(l)-Bi hexagonal layers, (d) Sketch of flux- 
grown crystals showing the axis [001]. (e) Recorded X-ray 
diffractograms. The upper trace (blue) is for powder speci¬ 
men of NasBi (1) crystals and the lower trace (orange) is for 
NasBi (2) single crystal. The inset in (e) shows a schematic of 
the Na-Bi phase diagram. The vertical dashed line indicates 
the Na-rich composition used in the flux growth. Adapted 
from Kushwaha et al. EE]. 


nodes, corresponding to the chiral anomaly sketched in 
Fig. [3j3 [IH3]t5HZl- As shown in Fig. [3]4, the lowest Lan¬ 
dau level (LL) of the Weyl states in a crystal is chiral. 
For x = 1 the lowest LL is right-moving (velocity v || B), 
whereas for x = — 1 it is left-moving. [IHH1I3E] At steady 
state, E transfers charge from say the x = — 1 branch to 
the x = 1 branch at the rate [5] 

The rate peaks for E || B and vanishes when E + B. 
In the weak-H limit, a Berry curvature approach gives 
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FIG. 5: Panel A: Sketch of the Landau levels (LL) in a 
Weyl semimetal showing chiral states in the lowest LL with 
opposite velocities and chiralities (arrows) || B. An E-field || 
B breaks chiral symmetry and leads to an axial current. Panel 
B shows the triangle anomaly that ruins the conservation of 
chiral charge. Panel C: The T dependence of the resistivity p 
in B — 0 and Hall coefficient Rh in NasBi. Rh is measured 
in B <2 T applied || c. At 3 K, Rh corresponds to a density 
n m 1.04 x 10 17 cm -3 . The inset shows the contact labels and 
the x and y axes fixed to the sample. Panel D: Curves of the 
longitudinal magnetoresistance p xx (B,T) at selected T from 
4.5 to 300 K measured with B || x and I applied to (1,4). 
The steep decrease in p xx (B : T) at 4.5 K reflects the onset 
of the axial current in the lowest LL. Adapted from Xiong et 
al. [27]. 


the chiral conductivity 0 


° x 47 i 2 hc c 


e 2 v ( eBv ) 2 


( 2 ) 


where r v is the intervalley life-time and ep the Fermi 
energy. 


IV. DIRAC SEMIMETAL NA3BI 

The Dirac semimetal NasBi grows as mm-sized, deep- 
purple, plate-like crystals with the largest face paral¬ 
lel to the a-b plane (c is normal to the planes) (Fig. 


[4^i, b). We annealed the crystals for 10 weeks before 
opening the growth tube. Details of the growth and 
characterization are reported in Kushwaha et al. [28]. 
To avoid oxidation, crystals were contacted using silver 
epoxy in an Argon glove box, and then immersed in para- 
tone in a capsule before rapid cooling. In NasBi, the 
Dirac nodes are located at the wave vectors (0,0, ±&d) 

with kp ~ 0.1 A [19j 20] . Initial experiments in our 
lab [26] on samples with a large Fermi energy ep (400 
mV) showed only a positive MR with the anomalous B- 
linear profile reported in CdsAs 2 [25] . 

Recent progress in lowering ep has resulted in sam¬ 
ples that display a non-metallic resistivity p vs. T profile 
and a low Hall density n ~ lx 10 17 cm -3 (Fig. HD 
We estimate the Fermi wavevector kp = 0.012 A (8x 
smaller than kp). The unusual profiles of p and the Hall 
coefficient Rp in Panel C imply the zero -B energy spec¬ 
trum shown in Panel B. Below ~10 K, the conductivity 
is largely due to electrons in the conduction band with 
electron mobility /i ^ 2,600 cm 2 /Vs). Because the energy 
gap is zero, holes in the valence band are copiously ex¬ 
cited even at low T. As T rises above 10 K, the increased 
hole population leads to a steep decrease in p and an in¬ 
version of the sign of Rp at 62 K. From the maximum 
in Rp at 105 K, we estimate that ep ~ 3kpT ~ 30 mV. 
As shown in Fig. HP, the resistivity p xx in a longitudinal 
field (B| |I, the current) displays a remarkable peak at 4.5 
K corresponding to a large negative MR (the resistance 
measured is R\a, 2 ?> {I applied to contacts 1 and 4, and 
voltage measured between contacts 2 and 3; the inset in 
Fig. |p shows the contact labels and the x and y- axes). 
Raising T above ^100 K suppresses the peak. The small 
density n implies that ep enters the lowest (N = 0) LL 
at B ~ 4-6 T. 


V. A NARROW CURRENT PLUME 

The axial current is predicted to be large when B is 
aligned with E. A valuable test then is the demonstra¬ 
tion that, if E is rotated by 90°, the negative magne¬ 
toresistance (MR) pattern rotates accordingly, i.e. the 
axial current maximum is selected by B and E, rather 
than being pinned to a crystal axis, even in the weak-T? 
regime. 

To test the anisotropy, we rotate B in the x-y plane 
while still monitoring the resistance ^ 14 , 23 - Figure [6jA 
shows the curves of the resistivity p xx vs. B measured 
at 4.5 K at selected 0 (the angle between B and x). The 
MR is positive for 0 = 90° (B| |y), displaying the nominal 

-linear form observed in CdsAs 2 [25] and NasBi [26] 
with B||c. As B is rotated towards x (0 decreased), 
the MR curves are pulled towards negative values. At 
alignment (0 = 0), the longitudinal MR is very large and 
fully negative (see SI for the unsymmetrized curves as 
well as results from a second sample). 

We then repeat the experiment in situ with I applied 
to the contacts (3, 5), so that E is rotated by 90° (the 
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FIG. 6: Evidence for axial current in NasBi (J4) obtained 
from transport measurements in an in-plane field B. Panel 
A shows plots of the resistivity p xx vs. B at selected field 
angles 4> to the a>axis (inferred from resistance #14,23, see 
inset). For 0 = 90°, p xx displays a B-linear positive MR. 
However, as 0 0° (B||x), p xx is strongly suppressed. Panel 

B shows plots of #35,26 with E rotated by 90° relative to 
Panel A (B makes angle (j)' relative to y; see inset). The 
resistance # 35,26 changes from a positive MR to negative as 
4>' —> 0°. In both configurations, the negative MR appears 
only when B is aligned with E. The spheres in the insets 
show the field orientations in 3D perspective. Adapted from 
Xiong et al. 


measured resistance is # 35 , 26 )- Remarkably, the observed 
MR pattern is also rotated by 90°, even when B < IT. 
Defining the angle of B relative to y as <//, we now find 
that the MR is fully negative when <f>' = 0. The curves 
in Panels A and B are nominally similar, except = 0 
and ft = 0 refer to x and y, respectively. 

We remark on the two striking features in Figs. and 
[6^3. The appearance of a strongly varying in-plane MR is 
rare in metals. The closest example we are aware of is the 
“anisotropic MR” (AMR) observed in ferromagnets (30]- 
[32] . In an in-plane B, one observes p(0) ^ P + Qsin 2 0 
caused by anisotropic scattering of carriers with velocity 
v||M vs. v _L M (M is the magnetization and P and Q 
are positive parameters). However, as AMR always leads 


FIG. 7: Magnetoresistance of NasBi (# 14 , 23 ) when B lies in 
the x-z plane at selected angles 0 , with E || x axis (6 is the 
angle between B and x). As in Fig. 6, an enhanced current 
(large negative MR) is observed when B approaches align¬ 
ment with E. When B is perpendicular to E, the MR has a 
positive H-linear profile. The field orientations are shown in 
3D perspective in the lower insert. From Jun Xiong (unpub¬ 
lished) . 


to a positive MR, it cannot account for the very large 
negative MR observed when </> (or (/)') 0. Moreover, no 

magnetic transitions have been detected in NasBi down 
to 2 K. 


VI. STEERING THE PLUME BY FIELD 

The second feature is the tight linkage between the MR 
pattern and the common direction of B and E even in 
weak B (<C 1 / p ~ 4 T). If the in-plane MR is caused 
by FS anisotropies (between a and b axes), the extrema 
of the oscillation ought to be pinned to the crystal axes. 
We should not be able to alter the resistivity tensor by 
rotating the weak E and B fields (this violates linear 
response). Viewed from ordinary transport, the observed 
rotation in weak B is indeed anomalous. However, it 
agrees with the prediction of the chiral anomaly - the 
axial current peaks when E aligns with B even for weak 
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FIG. 8: Angular dependence of the axial current in Sample 
J4 at 4.5 K inferred from measurements of i?i 4,23 in tilted 
B(0, <p). In Panel A, B lies in the x-y plane at an angle <p to 
x (sketch in inset). The conductance enhancement A a xx at 
fixed B is plotted against <p for fields 3 < B < 7 T. The insets 
show the polar representation of A a xx vs. <p. In Panel B, 
Aa xx is plotted versus 0 for 3 < B < 7 T for B lying in the x-z 
plane. As sketched in the inset, 0 is the angle between B and 
x. In Panels A (B), the axial current is peaked when <f> 0 

(0 —> 0) with an angular width that narrows significantly as 
B increases. Panels C and D show the orientations of B and 
E relative to the crystal axes in Panels A and B, respectively. 
Adapted from Xiong et al. El- 


fields. 

The large negative MR in Panels A and B suggests 
a long relaxation time for the novel current. We esti¬ 
mate the relaxation time r v for internode scattering from 
curves of the conductance G = l/i?35,26 as follows. At 
low T>, G increases rapidly as B 2 consistent with Eq. [ 2 ] 
We form the ratio cf x /(Jq = .(v/ T tr) (where 

a 0 is th e Drude conductivity, £3 is the magnetic length 
yJhfeB and r tr the usual transport lifetime). Fitting 
to the observed parabolic curve, we find that T v /r t r = 
40-60. The scattering rate relaxing the axial current is 


anomalously low compared with the scattering rate 1 / r tr 
of the conventional states in zero B. 

VII. ANGULAR WIDTH OF PLUME 


A surprise to us is the acute sensitivity of the novel 
current to misalignment at large B. We have exam¬ 
ined how the conductivity derived from ^ 14,23 decays 
as B is tilted away from x in either the x-y or the x- 
z plane. Figure [8^ displays the curves of A d xx (B,(/>) = 
a xx (B^(f)) — a xx (B, 90°) vs. <p as B is tilted in the x-y 
plane at an angle <p to x, with B fixed at values 3—>7 T. 
Figure [8^ shows the same measurements but now with 
B lying in the x-z plane at an angle 0 to x. In both 
cases, the low-field curves (B < 2 T) are reasonably de¬ 
scribed with cos p <p (or cos p 0) with p = 4 (not shown). 
However, for B > 2 T, the angular widths narrow signif¬ 
icantly. Hence, at large T>, the axial current is observed 
as a strongly collimated beam in the direction selected 
by B and E as <p or 0 is varied. 

To see what happens at larger F>, we extended mea¬ 
surements of Ri 4 ? 23 to H = 35 T. We observe a new fea¬ 
ture at Hk ~ 23 T when B||y. As B is tilted away from 
y {(j) 55°), the feature at becomes better resolved 

as a kink. The steep increase in p xx above suggests 
an electronic instability at large B. However, as we de¬ 
crease <p below 45°, Hj , c (</>) moves rapidly to above 35 T. 
The negative MR curve at p = 0 remains unaffected by 
the instability up to 35 T (the small rising background is 
from a weak B z due to a slight misalignment). 

To us, the unusual locking of the negative MR pattern 
in Figs. UK and Oh to E and B in weak B constitutes 
strong evidence for the axial current. The experiment 
confirms the B 2 behavior in Eq. [ 2 ] and provides a mea¬ 
surement of the long internode scattering lifetime. How¬ 
ever, the width of the collimated beam in the direction 
of B is much narrower than expected from the theory. 

In addition to the results here, several groups have also 
reported observing the chiral anomaly in other materials 
(primarily TaAs and ZrTes). The evidence is by and 
large restricted to the appearance of negative longitudi¬ 
nal MR often bracketed by large positive MR at lower 
and higher B. A concern is that a negative, longitudinal 
MR restricted to a narrow field interval is (by itself) a 
rather slender reed to hang a weighty claim from. Fur¬ 
ther tests, such as the demonstration of the field-steering 
effect, would appear to be necessary. Nonetheless, the 
dramatic increase in experimental activity on a growing 
list of candidate Weyl semimetals is an encouraging sign 
for the field. We anticipate exciting experimental devel¬ 
opments in the next few years. 
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